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Abstract 

We continue the study of the two-component charged Bose gas initiated by Dyson in 1967. 
He showed that the ground state energy for N particles is at least as negative as —CN'^/^ for 
large N and this power law was verified by a lower bound found by Conlon, Lieb and Yau 
in 1988. Dyson conjectured that the exact constant C was given by a mean-field minimization 
problem that used, as input, Foldy's calculation (using Bogolubov's 1947 formalism) for the one- 
component gas. Earlier we showed that Foldy's calculation is exact insofar as a lower bound 
of his form was obtained. In this paper we do the same thing for Dyson's conjecture. The 
two-component case is considerably more difficult because the gas is very non-homogeneous in 
its ground state. 
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1 Introduction 

In 1967 Dyson ^ showed that a system composed of non-relativistic, charged bosons is unstable 
in the sense that the ground state energy of N particles is at least as negative as -CN'^I^ instead 
of —CN, where C is some constant. A lower bound of the form —C'N'^^^ was derived later |CLYj . 
thereby establishing the correctness of the exponent 7/5, but not the constant C. 

In an earlier, parallel development, in 1961 Foldy F considered the problem of the one- 
component Bose gas ("jellium") in which charged particles (all of the same charge) move in a uni- 
formly charged, neutralizing background. Using Bogolubov's 1947 theory "B', Foldy "derived" the 
high density asymptotics for the ground state energy of this problem as proportional to —I^Np^l^ 
where /q is defined in Eq. (0) below. The correctness of this p^/^ law, but not the coefficient /q; 
was also proved in f jCLYj l. 

Dyson was motivated by Foldy's work, for he realized that if one treated one of the two com- 
ponents (say, the positive one) as a background for the other, and if one allowed the density to be 
variable, one would easily arrive — heuristically, at least — at the N"^/^ law. This will be explained 
below. 

Two obvious questions arise from this earlier work. What is the correct coefficient for the p^^^ 
law at high density and what is the correct coefficient for the A^^/^ law? The former question was 
resolved by us in |LSoj . where we showed that Foldy's Iq is, indeed, correct at high density as a 
lower bound. Foldy's calculation is "essentially" an upper bound, but some technical issues must 
be clarified. The proof that Foldy's calculation, indeed, gives an upper bound can be found in [S]. 

In ^ Dyson derives a rigorous upper bound for the N"^/^ law, but with a coefficient C that is 
clearly too small. He conjectures a "correct" coefficient, however, and in the present paper we shall 
show that Dyson's conjectured coefficient gives a correct lower bound (asymptotically as — > oo). 
An asymptotically correct upper bound for the two-component gas energy is also given in [Sj. 

Actually, our lower bound is slightly more general than just the case of A^/2 particles of each 
charge. We prove the lower bound for the case in which the total number is A^ without restriction 
to the A^/2, N/2 case. 

In order to understand the reason that the proof for the two-component case is more difficult 
than that for jellium, presented in |LSoj . it is necessary to recapitulate Dyson's argument briefly. 
His picture is that there is a local density of particles p(x), which has a local energy density given 
by Foldy's formula, i.e., Iqp{x)^^^. One might question whether the jellium energy can be simply 
taken over to the two-component situation, but it is correct to do so, as our lower bound shows. 
There is a good reason for this within the Bogolubov theory, as we shall explain in Sect.|3J but let 
us continue with Dyson's picture now. In addition to the local energy there is also a kinetic energy 
caused by the variation in p, namely / |Va/ p{x)\'^dx. Such an "envelope" energy is familiar from 
Thomas-Fermi- Weizsacker and Gross-Pitaevskii theories, for example. 

If this total energy is minimized with respect to p we are led to a differential equation in -^p 
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with the side condition that J p = N, but the basic features are clear. The scale length of p{x) will 
be of the order N~^/^, the amplitude of p will be N^/^ and the energy will be N^^. Indeed, if we 
define ^{x) = N'^/^ p{N-'^/^x) such that f ^"^ = 1 then the energy as a function of p that we 
have to minimize can be rewritten as 

U ivvpP-^o/ p'/' = Ny'Uf w-iJ ^^/^V (1) 

which makes the scaling explicit. We then have to minimize the right side with respect to <I> under 
the condition f = 1. The finiteness of this minimum energy is an easy consequence of the 
Sobolev inequality. The existence of a minimizing $ is a little harder and follows with the help of 
rearrangement inequalities. It satisfies a Lane-Emden differential equation for some /x > and all 
X 

- A$(x) - |/o $(x)^/2 ^ ^^>(x) = 0. (2) 

The uniqueness of <I>, up to translations, is harder still. See |Bel IkI ImSI IZ] . 

Dyson's conjecture is contained in |E1 eqs. 104, 105]. The normalization convention and units 
employed there are not completely clear, but the heuristics leading to (104) is clear. This will be 
clarified in Sect. 121 

The last topic to discuss in this introduction is the essential difficulty inherent in the two- 
component problem. As in the jellium proof in |LSoj we decompose space into suitable small boxes 
and impose Neumann boundary conditions on each. The scale length of these boxes is A^~^/^+^. 
Using the "sliding argument" of |CLYj we can ignore the Coulomb interaction between different 
boxes. We must then distribute the N particles in the boxes and what we must prove is that the 
number in the various boxes has a coarse grained density given by the solution p to Since the 
energy is super concave (i.e., —N"^^^), the lowest energy is obtained by putting all the particles in 
one box. What prevents this from happening is that the boxes are not really totally independent, 
thanks to the kinetic energy operator. In other words, we must somehow save a little bit of the 
kinetic energy operator to prevent wild variations in particle density between neighboring boxes. 

The conundrum is that the mean-field energy in uses all the kinetic energy, not just some of 
it. Likewise, to get the intrabox energy (the second term in (^) we also need the full kinetic energy. 
The resolution is to split the kinetic energy operator —A into a high-momentum part for use in 
calculating the intrabox energy and a low-momentum part for use in reproducing the first term in 
Naturally, error terms will arise and the chief difference between this paper and our earlier 
jellium paper is centered on the definition of the splitting and the management of the induced error 
terms. 

The proof of the main theorem starts in Sect. where we show how to localize the problem 
into large boxes of size L ^ N^^^^ with Dirichlet boundary conditions. This size is larger than 
the expected size of the bound complex. The decomposition of the kinetic energy into large and 
small momentum is carried out in an appendix. It is used in Sect. [51 which localizes further into 
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really small boxes of size i ^ N""^^^ . For the relevance of the scale N~'^/^ the reader is referred to 
the heuristic discussion in Sect. El The control of the electrostatics using sliding is also discussed 
in both of these sections. Sect. [HI discusses the ultraviolet and infrared cutoffs of the interaction 
potential with control of the errors they introduce. In contrast to the treatment in |LSoj we now 
need to use the kinetic energy to control the errors caused by the ultraviolet and infrared cutoffs in 
the potential. Sect. [7| controls all the unimportant parts of the localized Hamiltonian and reduces 
the problem to Bogolubov's Hamiltonian, which is analyzed in Sect. |SJ 

Sect, m gives the first, simple bound on kinetic energy, local non-neutrality, and an estimate 
on the local condensation. Sect. improves the estimate on condensation with the help of the 
method of "localizing large matrices" in |LSoj (and which is reviewed in an appendix). (Recall 
that in |LSoj we had to reach the final estimate on the various energies by a succession of finer 
and finer error bounds, each taking the previous bound as input.) In Sect. 1111 we give the final 
bound on the energy in each small box. We have to treat boxes with few particles as well as many 
particles separately. In Sect. we show how the kinetic energy estimate in the appendix (in which 
the low momentum kinetic energy between boxes leads to a difference energy on a lattice) leads, in 
turn, to the term J |Vy^p in the energy functional (^. In the final Sect. all the estimates are 
put together and we show how to choose the various parameters to get the desired minimization 
problem for the lower bound. 

We thank Y.Y. Li and M.I. Weinstein for pointing out the references |K1 IMSl IZ] to us. 

2 Basic definitions and main theorem 

We consider N particles with charges Cj = ±1, i = 1,...,N. The Hamiltonian describing the 
system is 

N 

//^=j:-iA.+ (3) 

i=l l<i<j<N ' 

acting on L^(M^^). We shall not specify the number of positive or negative particles, but simply 
consider the smallest possible energy 

E{N) = infjinf speci2(j^3N^H]\f \ a = ±1, i = 1, . . . , N}. 

Instead of considering Hj^ depending on the parameters ei,...,e7v we may consider it as one 
operator on the enlarged space L?' ((M^ x {1, —1})^), where the set {1, —1} contains the values of 
the charge variables. Then 

E{N) = infspeCi2((iR3x{i _i})iV)i?7V. 

Thus E{N) is the infimum of {'^ , H^"^) , over normalized functions ^' in ((M^ x {1,-1})^). 
Here we may restrict ^ to be non-negative and thus moreover also to be symmetric under the 
interchange of particles. Thus E{N) is the energy of a charged Bose gas. 
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Our main result in this paper is the asymptotic lower bound on E{N) conjectured by Dyson 
in The corresponding asymptotic upper bound is given in Together these results prove 
Dyson's formula. 

2.1 THEOREM (Dyson's formula). As N ^ oo we have 

E{N) = -AN'^/^ + o(iV^/^), 
where A is the positive constant determined by the variational principle 

0<«>, / $2<ll (4) 



with 



^ = inf<!i / \V1>\'-Io I <I>5/2 



It will be clear from the proof (see the discussion in Sect. that the error in Dyson's for- 
mula could have been written in the form N'^/^~^ for some ^ > 0. Although this is in principle 
straightforward we have not attempted to optimize the error term to determine the exact exponent 



3 Heuristic derivation of the energy 

In this section we give the heuristic derivation that leads to the local energy density I^p^^^ in the 
second term in (see also (0))). The fact that the constant for the two-component gas is the same 
as /q) the constant for the one-component gas, is somewhat mysterious. After all, there is not only 

the negative H — interaction energy, but there are also the -|--|- and energies. Moreover, the 

p in the two-component case refers to the total number of particles, which is twice the number of 
each charge (in the neutral case), whereas p refers only to the movable particles of one sign in the 
one-component case. 

Another problem is that the usual thermodynamic limit exists for the one-component gas only 
because we constrain the background to have uniform density and we do not allow it to contract to 
a high density ball, as it would if one merely minimizes the energy. The two-component gas does 
not have the usual thermodynamic limit because, as Dyson showed, its energy goes as 

Nevertheless, we shall describe the two-component charged gas on a small local scale £ as a gas 
of uniform density p. For such a uniform gas we shall imitate Foldy's calculation jF] to arrive at the 
energy —Iqp^/^1^ . We shall assume that £ ^ N~^/^ (the scale on which p{x) varies), but I ^ p^^/^ 
which, as we shall see, is the relevant scale for the uniform gas. Note that we expect p ^ N^l^ 
and hence ~ N^"^^^. In any event, we are being "heuristic" in this section and the reader is 

welcome to ignore this "derivation" if it is not pleasing. 
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We shall use periodic boundary conditions, as usual, and write the Hamiltonian in second 
quantized form in the manner of Bogolubov/Foldy as 

k 

1 \ ^ \ ^ 47r / ^ ^ * * 

+ 2£3 2^ 2^ "fc2 (,S+"g+"{g-fc) + "(p+fc)+ + %~"'q-'^{q-k)-"'(p+k)- 

ky^o pq 

-'^CL*j^a'*g„CL(g-k)~"'{p+k)+)^ (6) 

where the sums are over momenta in the set '^1? . Here a*_|_ creates a state with momentum p and 
charge ±1. 

In the second term we have excluded A; = 0, which is also what Foldy does. It would of course 
be meaningless to include this term and leave the Fourier transform of the potential equal to 
If, instead, we defined the Fourier transform of the potential for A; = to be proportional to 
(which is the maximal value of the Fourier transform for A; 7^ 0), the contribution from the term 
A; = would be proportional to —{p£^)i^^ for a neutral system. If, as we assume, i ^ p^^^^ 
then {pf)£''^ <C p^Z-^i^. Hence we may ignore the A: = contribution. In Foldy's situation the 
corresponding term would not contribute to the thermodynamic limit. 

The next step in the heuristic derivation is to exclude those terms in the second sum above 
that do not contain precisely two creation or annihilation operators of particles of momentum zero. 
Subsequently these zero momentum creation or annihilation operators Uq^. and Oq^ are replaced by 
the square root of half the particle number, namely p£^/2. We then arrive at 

Yl i^^ (4+afc+ + a-fc+«-fe+ + 4-afc- + a*^k-a-k-) 

k 

+ ^ ^ "p" + ^^-fc+«-fc+ + "fe+«-fc+ + "'k+(^-k+ 

k^O 

+ 'A-O'k- + O-fc-O-fc- + Ofc_alfe_ + ai^_a_j^_ 

'^l+O'-k- + ^*~k+^*k- + "fc+"-fc- + ^-k+^k-) ■ 

This expression should be compared with Lemma l<S.H where we arrive at a similar expression as a 
rigorous lower bound on part of the full Hamiltonian. (In comparing with the lemma one should 
replace 6^ a^^, iy^,iy~ p£^/2, Vr,R{k) 7^,4 1, (It^)^"^ 0, and of course also 

The final step is to recognize that the resulting quadratic Hamiltonian has the following property, 
as the reader can easily check: The operators a^_[_ always appear in the potential energy term (the 
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last sum above) in the combination = (a^_,_— a^_)/\/2- This is a normal mode since [d^, d*] = dk,q- 
The other normal mode = (a^^ + o-kJ)/ appears only in the kinetic energy term (the first sum 
above), i.e., the kinetic energy is a^_)_a^^ + al_aj^_ = d^Cf^ + d'^dj^. The ground state is achieved by 
having no c^, excitations, which leaves us just with the term d^d^ in the kinetic energy term. 

The conclusion is that the quadratic Hamiltonian is now exactly the same as Foldy's (but with 
df^ in place of a^) and, therefore, the ground state energy is —Iqp^^^I^. This conclusion could also 
have been arrived at by an explicit diagonalization of the total quadratic Hamiltonian. (See Sect.|SJ 
in particular Theorem 18.21 for comparison.) 

The detailed diagonalization analysis shows that the relevant momenta k are of magnitude p^^^ 
and hence as mentioned above the relevant length scale is The assumption £ ^ p^^/^ allows 

one to replace sums over the lattice by integrals J^i dk. This is how the integral in © 

appears. 

4 Localization 

If Dyson's conjecture is correct then the size of the boson cloud is proportional to A^~^/^. As a first 
step we shall localize the problem into cubes of size L, where we choose L as a function of N in 
such a way that N^/^L ^ oo as ^ cxD. Exactly how N'^/^L oo will be determined at the end 
of the analysis in Sect. 1131 As a consequence of our results we shall see that essentially all particles 
will concentrate within one of these cubes of size L. 

We shall do the localization in such a way that the cubes do not interact and the analysis can 
be done in each cube independently. The only thing to bear in mind is that the total number of 
particles in all cubes is A^. 

In analyzing the individual cubes we shall perform a further localization into smaller cells of a 
size i < L depending on A^ in such a way that N^/^i oo as N ^ oo (precisely how will again be 
determined in Sect. IT^ . 

We shall first describe how we control the electrostatic interaction between the different regions 
into which we localize. We do this in a manner very similar to what was done in |LSoj using the 
sliding technique of |CLYj . and shall use this technique both for the localization into the large cubes 
and again when we localize into the smaller cells. 

Let t, with < t < 1/2, be a parameter which will be chosen later in Sect. ^Jto depend on A^ 
in such a way that t — > as A^ — > oo. 

Let G,e G C7~(m3) satisfy 

1. < 6,0 < 1, e{x) = e{-x), and e(x) = e(-x). 

2. supp C [(-1 + t)/2, (1 - ^)/2]^ supp e C [(-1 - t)/2, (1 + t)/2f. 

3. 9{x) = 1 for X G [(-1 + 2t)/2, (1 - 2^)/2]^ and e(x) = 1 for x G [(-1 + t)/2, (1 - ^)/2]^ 
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4. All derivatives of order m for m < 3 of the functions 6, VT^-W, Q are uniformly bounded by 
Ct~^, where C is some universal constant. 

5. For all x G j^^ve ©(^ - = 1- 

We introduce the two constants 7, 7 such that 7 J 6{y)'^dy = 1 and 7 J @'^{y)'^dy = 1. Then 

l<7<(l-2^)"^ + <7< (1-*)"^- (7) 

We also introduce the Yukawa potential Ym{x) = |x|~^e~™l^l for m > 0. For m = this is, of 
course, the Coulomb potential. 

4.1 LEMMA (Electrostatic decoupling of boxes using sliding). 

There exists a function of the form uj(t) = Ct~'^ (we assume that u(t) > 1 for t < 1/2) such that 
for all xi,X2, . . . ,xn G IR'^, all ei, 62, . . . , cat, with |ej| = 1 for i = 1,2, . . . , N , all m > 0, and all 
X > we have 

eiejYm{xi-Xj)>-i I ^ eiCjO - Y^^^{x,-Xj)e -z^dz- ^'^^^^ (8) 
l<i<i<Ar ^3 l<i<j<N 

and likewise 

l<i<j<N ^3 l<i<j<N 



^ eiejYm{xi -Xj)>j ^ eje^G^ (^y - ^^.^^^(a;* - Xj)Q^ - dz 



2X 

(9) 



Proof. Since 9 and 0^ have the same properties it is enough to consider 9. By rescaling we may 
assume that A = 1. We have that 

j9{x - z)Ym+u;{x - y)9{y - z) dz = h{x - y)Ym+u;{x - y), 

where we have set h = j9*9. We chose 7 such that 1 = /i(0) =7 / 9{y)^ dy. Then h satisfies all the 
assumptions in Lemma 2.1 in |CLYj . We then conclude from Lemma 2.1 in |CLYj that the Fourier 
transform of the function F{x) = Ym{x) — h{x)Ym+oj{x) is non-negative, if to is chosen large enough 
depending on h. [The detailed bounds from |CLYj show that we may in fact choose to = Ct~'^, 
since uj has to control the 4th derivative of h.] Note, moreover, that lim^^^o -^(2;) = u^- Hence 

^2 AT 

e^ejFiy, -y^)>-^^ = _. 

^<i<j<N i=l 

The lemma follows by writing Ym^x) = F{x) + h{x)Ym+u){x). □ 

We shall use the electrostatic decoupling ® when we localize into the large cubes of size L and 
then Q when we localize further into the smaller cells. We begin with the localization into large 
cubes. 
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4.2 THEOREM (Localization into a large cube). Let 



N 



Hn,l = ^ + ejej-jY^jt) {xj - xj) 



l<i<j<N ^ 



be a Hamiltonian acting in the space 

L2 (([-L/2, L/2]3 X {1, -1})^) . Here Ad refers to the Dirich- 
let Laplacian in the cube [—L/2,L/2]^. Let 

El{N) = inf speci^AT,!,- (10) 

Then 

E{N) > El{N) - (ct-'^{N^'^Ly^ + \uj{t)N~^'^{N^/^L)-^^ . 

Proof. For all z G let ©z(x) = Q{{x/L) — z). We consider 0^ as a multiplication operator on 
L2(m3). We have for all z e that 



Using (jUJ with m = and that Jj^g f{z)dz = J2qez^ -^[-1/2 1/2]=^ -^("^ ~'~ '^^'^^ ^^^^ 

E 



> X] / ^k,2(2;i, . . . ,XAr)2 E h,k,eiej^Y^(t)iL{xi- Xj)dz 

where 

Fk,2(xi, . . .,xn) = QkT,+z{xi) ■ ■ ■ QkN+zixN), k = (fci, . . . , /cat) E Z^^. 
We localize the kinetic energy using the formula 

-A=Y^ e,+,(-A)e,+, - E (VG,+,)2 > e,+,(-A)G,+, - c(tL)-2, 

which gives 

N N 

^-A, > E / Fk,^(xi,...,Xiv)E(-Ai)FM(xi,...,Xiv)d2-CAf(tL)-2. 

i=l kGZ3iVj_^/2^^/2]3 ^=1 

The previous estimates imply that 
Hn> Yl I i^k,.(xi,...,X7v) 5^^£i^k,.(xi,...,X7v)dz-^^-CiV(tL)-^ 

k6Z^^[-l/2 l/2]3 9eZ3 

where we have introduced the operator 

AT 

^£ = Yl -I'^fc.-g^!? + Yl ^k,,q^k„qeiejjY^^t)/L{xi - xj) 

1=1 l<i<j<N 
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where Aj^"^^-* denotes the Dirichlet Laplacian in the cube {{q + z)L} + [— L/2,L/2]^. Note that 
the operator H^^^ above acts on functions for which the space variables are in the set {{q + z)L} + 
[— L/2,L/2]^. The operator H^^], is hence unitary equivalent to H^^l with replaced by 

N,{k) = #Mgik), M,{h, ...,kN) = {i\ ki = q}. 

We therefore get the lower bound 

Hn> [ ^m(xi,...,x,v)' j;i^L(iV,(k))dz-^^-C7iV(tL)-2. 

The theorem follows if we can prove that for fixed k E 

EL{N,{k)) > El{N). (11) 

To conclude (|lll) . let tpg be a normalized A'^g-particle ground state eigenfunction of -^at^.l, such 
that 

(A normalized 0-particle function we consider simply to be the number 1.) Define the A'^-particle 
eigenfunction 

<? 

(for fixed k there are only finitely many factors different from 1 in the above product) where we 
shall choose the (finite) parameters Tg £ {1,-1}. We then find that 

geZ3 q^q' 

(note that that there are again only finitely many non-zero terms in the sums) where 

Iqq' = (i^qi^q',Y'^^''^'1^''j'l'^^^i^^'^(^)/I^^^^ -^Mq'^q 
i j 

It is now easy to see that by an appropriate choice of the finite number of parameters Tq, we can 
make sure that the last sum in (fT^ is non-positive. In fact, if we average over all possible choices 
of = ±1 the last sum in ((T^ averages to 0. 

The estimate (|11|) is now a consequence of the variational principle. □ 

The rest of the analysis is concerned with estimating the energy El {N) from (|10|) . As explained 
in the beginning we may define El{N) as inf speci^iv^L when we consider -^at^l as an operator on 

N 

the symmetric tensor product <^gL'^{[—L/2, x {1,-1}). 

N 

In the space (^^L^(M'^ x {1, —1}) we shall use the notation of second quantization. If u G 

oo N 

L^(M^ X {1,-1}) then a*{u) denotes the creation operator in the Fock space (^^L^(M^ x 

Ar=0 
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{1, —1}). Products of the form aQ{u)aQ{u') or ag(M')ao(u) are however bounded operators on each 

N 

of the fixed particle number spaces x {1, —1}). We may consider Hj\f^L as acting in this 

space but restricted to functions with support in ([— L/2, x {1,-1})^. 

If n G we use the notation aj_(n) for the creation operator which creates the function u 

with charge ±1 respectively, i.e., the function u defined on x {1} or x {—1} respectively. 

5 The localization of the operator Hjy ^ 

We turn to the second localization into smaller cubes of size i > 0, where we shall choose £ as a 
function of N in such a way that £ < L and N^/^l oo as N ^ oo. As mentioned the precise 
choice of ^ as a function of N will be made at the end of section when we complete the proof of 
Dyson's formula. 

Localizing the kinetic energy is difficult. In the paper |LSoj we used Neumann boundary condi- 
tions on the small cubes. In this way we got a lower bound by studying independent small boxes. 
In the present situation we will not get the correct answer if we bound the kinetic energy below 
by Neumann energies in independent boxes. As explained the "kinetic energy" between boxes is 
important. A lower bound on the kinetic energy which both contains a contribution from within 
each box and a contribution that connects boxes is given in Appendix lEl We shall use this result 
now. 

We use the function 9 to do the localization. For z G we define Xz{x) = 0{{x/£) — z). Then 
supp Xz C {zi} + [(-1 +t)^/2,(l -t)£/2]^ Let r?^ = ^/T^xf. Then ||a"X^||oo < C(^t)-I"l and 
||<9°%||oo ^ C(£t)^l°l for all multi-indices a with \a\ < 3. 

Let Xz be the characteristic function of the cube {zf\ + [— ^/2, ^/2]'^, i.e., Xz{x) = Xq{x — iz). 
Let Vz denote the projection onto the subspace of {{{zf\ + [—1/2,1/2]^) x {1,-1}) consisting 
of functions orthogonal to constants. We shall consider Vz as a projection in L^i^ x {1,-1}). We 
define for each z G the operator 

The operators /C^^-* shall play the role of the kinetic energy within each box. 
Let aQj^{z) be the creation operators 

al^{z)=al{r'"'Xz). (14) 

i.e., the operator creating the constants in the cube {zf\ + [—1/2,1/2]^ . We introduce the notation 

Tv^ = aQj^{z)aQj^{z) and nz = n^ + n'^. (15) 

We refer to as the number of particles in the condensate in the cube {zl^ + [—1/2,1/2]'^ . The 
operator for the total number of particles in the cube {z€\ + [—1/2,1/2^ is 

N 

= ^X,(x,) (16) 

i=l 
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The operators for the total numbers of positively and negatively charged particles are determined 
by 

N 

+ K = '^z and -v^ = ^^ eiX^{x,i). 

i=l 

We refer to z?^ — as the number of excited particles. 

The kinetic energy that connects boxes will be a type of lattice Laplacian. In fact, for a map 
S : ^ M we define a lattice Laplacian 

^(^)= E YA^S{a^)-S{a2)f. (17) 

|(T]^ — (72 I — lo-^— (72| — 

This specific form of T has been chosen because it is convenient when we shall later compare with 
a continuum Laplacian (see Sect. [T^. 

Using Theorem inU in Appendix [0 with Q = [-{L + e)/{2£),{L + £)/{2e)]^ we obtain the 
following result. 

5.1 LEMMA (Kinetic energy localization). If ^ Li^iM? x {1,-1}) is normalized and 
has support in ([-L/2, L/lf" x {1, -1})^ then for all e > 0, < t < 1/2 

{i+e+ct')U,f2-A,^] > [ UY.f:{^t''^-'^^:^)A+ns!)dz 

\ ^=l / ,g[_l/2,l/2]3 V -&'^=^ J 

where — A^^^ is the Neumann Laplacian for the cube {z£} + [—£/2,i/2]^ and Sf : ^ M zs the 
map 

sficT) = [^J{nt+^ + n-+.) + 1 - 1) = (V(^.+<x) + 1 - 1) , (18) 

where 

In order to arrive at this lemma from Theorem IB . 1 1 we have used the inequality 

T{sT)<{T{st) + ns;)), 



where Sfia) = £-'^ ^nf^^ + 1/2. To prove this inequality first note that by the 2-dimensional 
triangle inequality 



'Sf + Sf - ^Sl + sij < {Si - S2f + {Si - S2f, (19) 

for all Si, S2, Si, S2 G IR. The estimate (|19p implies, in particular, that 

n^T{e~^Vn + 1-1) 
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is a convex map on non-negative functions n : ^ M. Hence 

T{sf) = TirW{nz+a) + i- 1) < (r(r V^.+. + 1 - 1)) < {T{st) +t{s-)), 

where the last inequahty is again a consequence of ((T^. 

As before we use the shding technique to locahze the electrostatic energy. If we combine 
Lemma l4. II and Lemma l5. II we obtain for ^' S x {1, —1}) normalized and with support 

in ([-L/2,L/2]3 x {1,-1})^ that 




2G[-l/2,l/2]3 

-CL^/£^ - uj{t)N ( ^ - ^ ) (20) 
where 

7e,t = (l + e + Ct3)-i(^7)-i (21) 
(recall that 7 and 7 also depend on t) and for z G we have introduced the operator 



N 



H. = E[-f'^i-^^^rieu)+ E e,e,w.(^,,x,), (22) 

1=1 l<i<i<N 

with 

Wz{xi,Xj) = Xz{xi)Y^(t) {xj -Xj)Xzixj). (23) 

e 

We shall choose e depending on N in Sect. in such a way that e ^ as ^ 00. We shall 
in particular assume that e < 1. Hence we may assume that ^e^t is bounded above and below by 
constants. 

The expressions 

Eaezs +l-fe,tT{Sf) are equivalent for the different z G [-1/2, 1/2]^. 

It is therefore enough to consider z = 0. Moreover, the operator X^o-gz^ -^o- commutes with the num- 
ber operators z/^ (the operators giving the number of positively and negatively charged particles 
in the cube {ia} + [-i/2,e/2]^). Note that if ^ has support in the set {[-L/2,L/2f x {1, -1})^ 
then for z outside the set f^ = [~^~^)^ + we have that ^' is in the kernel of the operators 

The estimate H20p then implies the following result. 
5.2 THEOREM (Localization into small cubes). With the definitions from 17 ^) . 173^7^) . and 



we have that 

El{N) > 77^ inf^^ 1^ 1^^- E + ^^(-^o*) j - CL'/e' - u;{t)N 



1 1 

2£~ 2L 



where the infimum is over normalized functions 

N L I L I 3 

G Wo = G <S)s^^(^' ^ [-1' M ^^"^ = ° ^^'^ ^ ^ ^ ["2£ " 2' 2£ + 2] }' 
and < t < 1/2 is the parameter from the beginning of Sect. ^ 
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In the following Sects. I ^lll we shall study the operators H^j, a G 'L^ . Since H^j commutes with 
we may simply restrict to the eigenspaces of z?^. We shall therefore not think of and v^r 
as operators, but as classical parameters and Vfj. Which values of these parameters that will 
give the optimal energy is of course not known a-priori, but they must satisfy v^j = N. The 
operator H^j is in written as an A^-particle operator, but it depends only on z^o- particles. We 
shall therefore simply think of it as a fo-particle operator. The operators H„ for different a are all 
unitarily equivalent. We shall in Sects. I^HTTl simply omit the subscript a. 



6 Long and short distance cutoffs in the potential 

Our aim in this section is to replace w in (|23() (omitting the index z) by a function that has long 
and short distance cutoffs. 

We shall replace the function w by 

WrM^^ y) = X{x)Vr,R{x - y)x{y) (24) 

where 

_ -\x\/r 

Vr,R{x) = Yr-. {x) - y,-l (X) = (25) 

Here < r < < ijj{t)^'^t. Note that for x < r then Vr^R{x) r^^ - R^^ and for \x\ » R then 
Vr,R{x) ^ |x|-ie-l^'l/^. 

We first introduce the cutoff R alone, i.e., we bound the effect of replacing w by WR{x,y) = 
X{x)Vr{x — y)x{y), where Vr{x) = |x|~^e~l^l/^ = Yji-i{x). Thus, since R < uj{t)^^i, the Fourier 
transforms satisfy 

(We use the convention that f{k) = f f{x)e~^''^dx.) Hence 

w{x, y) - wr{x, y) = x{x) {Y^/e, - Vr) {x - y)x{y) 

defines a positive semi-definite kernel. Note, moreover, that {Y^j^ — Vr) (0) = i?^^ — uj/l < R~^. 
Thus, 

eiej{w{xi,Xj) - WR{xi,Xj)) = ^ eiej{w - WR){xi,Xj) 

l<i<j<N l<i,j<N 

N 2 

-"^^{W - WR){Xi,Xi) 
i=l 

> -\v{Y^/i-VR){S))>-\vR-\ (26) 

Here we have used that '}2d=i X{xi)^ < J2iLi -^{^i) = ^• 

We shall now bound the effect of replacing wr by Wr^R- I.e., we are replacing Vr{x) = 
|2,|-ig-|a:|/_R |-jy (g-^'IZ-R _ e"'^'/*"). The correction is 1^-1(3;) = |x|~"^e~l^l/''. 
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6.1 LEMMA. We have for all 6 > the operator inequality 



N 



IT 

+ E e,ejX{xi)Yr~i{xi - xj)x{x,) > -Cv+u-{r''r^ + 5-^'\^'^). 

i=l l<i<j<N 

Proof. We set z = and omit the index z in this proof. For D > we define 



Wd{x) 



Y^-i{D), \x\<D 
y^-i(x), |x| > D 



Then the measure 

fiD = {47T)-\-A + r-^)WD = (4vr)-V-2y,,-i(Z?)l|,|<^- (47r)-iy;_i(Z?)5(|x| -Z^) 

is non-negative and Wd = fJ-D * Here refers to the radial derivative of l^-i and l|2:|<z) 

is the characteristic function of the set {|x| < D}. 
Define the (signed) measure 

TV 

d/x(x) = ^ eiX{xi)dfiDi (x - Xi), 

i=l 

where 

A = lmm{\xi-Xj\ \j = l,...,N, aej = -1, Xi € [-i/2,l/2f} , 

is half the distance from the particle at Xi to the nearest particle of opposite charge in the cube 
[—i/2,£/2]'^, i.e., the cube in which X is supported. Thus for Xi,Xj G [—£/2,£/2]^ we have 



Y,-i{x - y)diJ,D,{x - Xi)dfiDj{y - Xj) = Wo^ix - Xj)diJ,D,{x - Xi) 



< J Yr~i{x - Xj)dflDi{x - Xi) 

= Woiixi - Xj) < Y,.^i{xi - Xj) 
and both inequalities become equalities if CiCj = — 1. Hence 

eiejX{xi)Yr~i{xi - Xj)x{xj) > \ j j Y^.-i{x - y)d^i{x)d^i{y) 

l<i<j<N 



k^^i^i)'^ / / Yr-i{x - y)dfiD,{x)dfJ-DAy) 
i=l 

N 

> -l^Xixif / / Y^^i{x - y)diiD,{x)dnDAy) 



since Y^.-i is positive type. This inequality is very similar to Onsager's electrostatic inequality [Uj . 
We calculate 

j I Y,-i{x - y)di^D{x)d^iD{y) = J WD{x)dtiD{x) = lr-^D^Y,.-i{D)^ - D^Y;.,{D)Y,-.{D) 
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We have therefore proved the operator estimate 

^5 ^5 

-^^i.Ncu + ^ eiejX{xi)Yr-i{xi - Xj)x{xj) > ^ --Aj^Neu - Vi, 
i=l l<i<j<N i=l 

where for i = 1, 2, . . . , A'^ 

Viixi) = -Ixixif {D,r-^ + 3i?ri ^ 3^-1^ ^^2D./r _ 

Note that here Di depends on Xi and the positions of all the particles with opposite charge of Cj. 

From the Sobolev estimate for the Neumann Laplacian in a cube of size ^ we have the general 
lower bound 

-ANeu -V>-C j _ (^£-3 j y 

This gives for = ±1 that 

2 J]R3 \ \x\ rj \r'^ \x\ r J 

> _c7i.T(r3/VV2 + rV). 

When summed over i this gives the result of the lemma. □ 
If we combine the bound (|26|) and Lemma IH. II we have the following result. 

6.2 LEMMA (Long and short distance potential cutoffs). 

For < 5 < e < 1 consider the Hamiltonian 

N 

Hr,R = E(^^f^-^(^-^A^Seu)+ E e.e,WrA^,,x,) (27) 

i=l ^<i<j<N 



where Wr,R is given in \24\j and \25\) . If < r < R < Lo{t) ^£ then the Hamiltonian H defined in 
l^22\) obeys the lower bound 

H > -\vR-^- Cv+u- (r 3/2r + r V) . (28) 

IfO<r<R = uj{t)-^l we get 

H > - Cu+v-{5-^l^r^l^ + rV2). (29) 

A similar argument gives the following result. 

6.3 LEMMA. With the same notation as above we have for 0<r<r'<R'<R< uj{t)^^i that 

Proof. Simply note that 14', R' {x) — Vr^R{x) = 1/?/-! (x) — Yr-i (x) + Y^-i (x) — Y^i~i (x) and now use 
the same arguments as before. □ 
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We shall not fix the cutoffs r and R nor the parameter 6, but rather choose them differently 
at different stages in the later arguments. Moreover, we shall choose r and R to depend on v. We 
give an example of this in the following corollary. 

6.4 COROLLARY (Simple bound on the number of excited particles). 

For any state such that the expectation {H) < 0, the expectation of the number of excited particles 
satisfies {v — n) < i' m.m{l, Ce~^ ii^^^^} . (Recall that H depends on e.) 

Proof. We simply choose r = R and 6 = e/2 in Lemma 16.21 Then > C{e/2){u — n)i~'^ since 
{u, —A^l^u) > 7r^£~^||n|p for functions u in the cube {z£} + [—1/2,1/2]^ orthogonal to constants. 
We therefore have that 

H > C{e/2){u - n)r^ - \ur-^ - C7i/+i/- (£-3/2^1/2 ^ ^-3^2^^ 

Strictly speaking Lemma 16.21 requires r < R < uj{t)^^£. If however R = r > uj{t)~^£ we would get 
an even better estimate than the one stated if we set r = R = (jj{t)~^£ and 6 = e/2. Since in this 
case we may ignore the error — ^z/r^^. 

We now estimate u^v" < and make the explicit choice r = min{ez^~^/^, ^z/^"*^/^}. When 

eV 

-2/3 < ^j,-i/3 

I.e., V 

1/3 > ££-1 

we obtain 

H > C{e/2){iy - n)£''^ - Ce^^v^''^ - Ce^£-^v^l^ > C{e/2){v - n)£'^ - Ce'^v^'^ 
and when ez^~2/3 > £jy-i/3^ ^ jyi/3 ^ ^^-i obtain 

H > C{e/2){v - n)£-^ - C7e-3/2^i/2^ii/6 _ ^7^-^4/3 > C{e/2){v - n)£-^ - CrV^/^ 
These two bounds give that 

{v-n) < CvTn<aL{{e~^£v^/^f,e-^£v^/^}. 
We of course clearly have that v — n < v. We therefore get the result claimed. □ 

7 Bound on the unimportant part of the Hamiltonian 

In this section we shall bound the Hamiltonian ^ given in (|27|) . We emphasize that we do not 
necessarily have neutrality in the cube, i.e., and v" may be different. We are simply looking for 
a lower bound to that holds for all v^. The goal is to find a lower bound that will allow us to 

conclude that the optimal value of — v~ , i.e., the value for which the energy of the Hamiltonian 
is smallest, is indeed close to zero. We shall also conclude that the number of excited particles 
— n is small. These conclusions will be made in Sects. l^and lTUl 
In this section we shall fix ^ = 0, i.e., we are working in the cube [—£/2,£/2]^. We shall simply 
omit the index z. 
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We shall express the Hamiltonian in second quantized language. This is purely for convenience. 
We stress that we are not in anyway changing the model by doing this and the treatment is entirely 
rigorous and could have been done without the use of second quantization. 

Let Up, ip/ir G (N U {0})'^ be an orthonormal basis of eigenfunction of the Neumann Laplacian 
— Anbu in the cube [—£/2,i/2]^. More precisely, — Anbu^^p = IPp'Up- I-e., 

3 

Up{xi,X2,X3) = Cp£"^/^ cos {Pj{xj + i/2)) , 

where the normalization satisfies cq = 1 and in general 1 < Cp < \/8- The function no = ^-3/2 is 
the constant eigenfunction with eigenvalue 0. We note that for p 7^ we have 

(up, -Ancu^^p) > vr^^"^. (30) 

We now express the Hamiltonian from (|27() (omitting z) in terms of the creation and 

annihilation operators a*^ = a±{up)* and Cp^ = a±{up) . 
Define 

^pg.Mi/ = y" j Wrfi{x,y)up{x)uq{y)u^{x)uy{y) dxdy. 
We may then express the two-body potential in -ff^^ as 

Motivated by Foldy's and Dyson's use of the Bogolubov approximation it is our goal to reduce 
the Hamiltonian so that it has only what we call quadratic terms, i.e., terms which contain 
precisely two with p 0. More precisely, we want to be able to ignore all terms in the 
two-body-potential containing the coefficients 

• ""^00,00- 

• ^pO,qO ~ ^Op.Oq' where p,q ^ 0. These terms are in fact quadratic, but do not appear in the 
Foldy Hamiltonian. We shall prove that they can also be ignored. 

• ^po.oo = ^Op,oo = ^^oo.po = ^oo,Op> where p^O. 

• ^P9,/.o = %o,pg = ^qp,o^, = ^0/.,gp' ^here p,q,n^ 0. 

• Wpg^f,^, where p, q,^i,u 0. 

We shall consider these cases one at a time. 

7.1 LEMMA (Control of terms vi^ith i«oo,oo)- 

The sum of the terms in containing Wqq qq ^-^ equal to 

^■"^oo.oo [(^^ - - n+ - n-] . 
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Proof. The terms containing u;qq qq are 

^'"^00,00 (^0+^*0+ '*o+'^o+ + '^o- ^0-^0-^0- ~ 2ao+ao_ao_ao+) , 
which gives the above result when using the commutation relation [a^, a*] = 6p^q. □ 
7.2 LEMMA (Control of terms vi^ith WpQ^go). 

The sum of the terms in ji containing WpQ^Q or w^pQ^ with p,q ^ is bounded below by 

-4TrR'^r^\n+ - n~|(z/ - n). 
Proof. The terms containing WpQ^gQ = w^p Q^ are 

^|c^ \ 



= 2^^pO,qo("^ - n - «p-S 

Note that commutes with ^ '"^po,qO*^p±^q±- 

D 

^ ^'^ 1 1 dyup{x)uq{x) dx. 



p^O 

q^o 

We have that 



■"^pO,qO 



Hence 



X] ^pO,goap±aq± ^ ^ ^ Wr,R{x,y)dy ^Up{x)a*p^ X?ip(x)a*± dx. 

go Vp^o y \p^o J 

< £~^sup / Wr,R{x',y)dy / y^ Up{x)al^ y^ Up{x)at^ dx. 

= £~^sup / Wr,R{x',y)dy^alj^ap^. 
^' p^O 

Since X^^^q "pi^^pi = i^"*" - and 

sup y Wr,R{x,y)dy < J Vr,R{y) dy < 4ttR^ 
we obtain the operator inequality 

p^O 
q^O 

and the lemma follows. □ 

Before treating the last three types of terms we shall need the following result on the structure 
of the coefficients w„„ 



20 



7.3 LEMMA. For all p', q' E {tt/£) (N U {0})^ and a G N there exists Jp^, E R with J°,^, = J^^^, 
such that for all p, G (^/^) (I^ U {0})^ we have 



JpiiJqu- (31) 

Moreover we have the operator inequalities 

0< ^ %p',ooap±V±= Wpo^Qp,a*p±ap,^<4TTr^R^{iy^-n^) (32) 

and 

< J]] Wp^,mp'a*p±ap,± < r-^{v^ - n^). 

p,p',m^O 

Proof. The operator A with integral kernel Wr,R{x, y) is a non-negative Hilbert-Schmidt operator on 
L2(m3) with norm less than sup^, K-,i?(fe) < 47ri?^. Denote the eigenvalues of ^ by A^, a = 1, 2, . . . 
and corresponding orthonormal eigenfunctions by ipa- We may assume that these functions are 
real. The eigenvalues satisfy < < AttR"^. We then have 

= ^A„ / Up{x)Uf,ix)ipaix) dx / Uq{y)u^{y)ipo,iy)dy. 

1 11 

The identity pT|) thus follows with Jp^ = Xj J Up{x)u^{x)ipa{x) dx. 

If as before V denotes the projection orthogonal to constants we may also consider the operator 
VAV. Denote its eigenvalues and eigenfunctions by and ip'^. Then again < A^ < AttB?. Hence 
we may write 

^po.op' =^'^Y^'a j Up{x)ip'^{x) dx / Up>{y)ip'^{y)dy. 
Thus since all are orthogonal to constants we have 

p,pVo 



^"'E^'" (E / n,{x)^'^{x)dxal\ W [ 



Up(x)9J^(x) a. 



'p± 



The inequalities (|,32j) follow immediately from this. 

The fact that Yip p' mj^o''^pm mp'^p±^p'± — *-* follows from the representation (|^ . Moreover, 
since the kernel Wr,R{x,y) is a continuous function we have that Wr,R{x,x) = X^a -^aV'a(^)^ 
almost all x and hence 
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We therefore have 



< suY)Wr,R{x' ,x') / ^^tp(a:)a*± ^^ip(a;)a*± dx = svL\)Wr,R{x' ,x'){u^ -rr^) 

Vpv^O / \P7^0 / 



and the lemma follows since sup^./ Wr,R{x' ,x') <r. □ 
7.4 LEMMA (Control of terms with Wpg^^^). 

The sum of the terms in ^ containing Wpq^^y, with p,q, fj,,iy is bounded below by 

1^{Y1 Jp^,{a*p+a^^ - a*p_^a^_)J - ^r'^u - n). 

Proof. The relevant terms are 

pq,^iv^O 

= h ^ ^pq,tiu {a*p+af,+al^a^^ + a*p_a^_al_a^_ - 2a*+a*_a^_a^+) 

pq,IJ.uj^O 

_1 ^ * _ 1 ^ * 

2 / > ^pm,mp''^p+^p'+ 2 / j ^pin,mp''^p~'^p' ~ 
pp' ,m^0 pp',m^O 

2 



^5Z(Z] '^pmK+"m+ -"p-^A'-)) -I'^ \i^-n} 

a p/^7^0 



where we have used the last estimate in Lemma 17.31 and that z^"*" — + u —n = u — n. □ 
7.5 LEMMA (Control of terms vi^ith u?po,oo)- 

The sum of the terms in ^ containing w^^ qq, ^'op,oO''"^oo,pO' ^oo,Op' p is, for all 
e' > 0, bounded below by 

- £'-^A-Kr^B?{n + - n) - 2e'wQQQQ ((n+ - n'f + l) , (33) 

or by 

(z^+ - V-) ^ u)po 00 (op+ao+ - o*_ao- + oS+ap+ - 00-^^.) 
-e'-Mvrr^iJ^^^ + _ fi) _ 2e'woo 00 ((1/ - n)^ + l)). (34) 
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Proof. The terms containing u)p0,00' ^Op,oO'^oo,pO' ^oofip are 

— 2ag_ao_ao+%+ ~ 2ao_ao-^p+^o+ ~ ^ag+ao+^o-^p- ~ 2ao^ao+ap_ao_^ 
= ^ ^ J^O'^oo ((""^ - ^")(a;+ao+ - ap-«o-) + («o+«p+ ~ ao~ap-)(^^ ~ ^")) 

a p^O 

- ^ ^ J^o-^oo (a;_ao-(^^ - n" - 1) + a*o_ap_(n+ - - 1)) 

a p^O 

> E (^oo)' ((^^ + If + (n+ - n- - l)^) 

a 

Here we have used the representation and in the last step a Cauchy-Schwarz inequahty. We 
arrive at the first bound in the lemma since from (1321) we have that 



JpoJp'o(h>±(^p'±= Y ^po.op' ap±ap'± < 47r£ ^R^{v^ -n^). 

a p,p'^0 P,pV 

The second bound l\M\i follows in the same way if we notice that the terms containing t5po,oO' 
''^Op,oO'^oo,pO> or Woo.Op ™ay be written as 

+ E E J°o Joo ( (^^ -i^^ -n- + iy-){a*p^ao+ - al_aQ_) 



p^O 



□ 



7.6 LEMMA (Control of terms with Wpg^rno)- 

The sum of the terms in H^^^ containing iVp^^^Q, %p,om'^mO,pg' ^'^ ^Om,qp' "^^^^ p,q,m ^ is 
bounded below by 

2 



for all e' > 0. 

Proof. The terms containingWp^ , w^p^^^w^Qp^, or w^^^p, with p,q,mj^O are 

^pq,mO ( ^p+^g+'^m+O'0+ + O'0+a'm+O'g+ap+ + ap-Oq-am-'^O- + ^O-'^m-Og-O'p- 

pqm^O 
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"^p+^g-'^rra+^O- ^0-^m+'^g-'^p+ "p-^g+^m-^0+ %+'^m~^q+%- 



a (j^O pmj^O 

+ Jpm{^p+^m+ ~ ^p-^m-) •^go(^g+'^0+ ~ ^g-'^O-) 

o pm^O q=iO 

I ^ ^ ^ ^ ^ / ]^ ^ ^ ^ ^ 

> — e ^ ^50,0g'^5+'^q'+'^0+'^0+ ~ ^ 2-„ '^gCOq'^g-'^q'-^O-'^O- 

2 



-2e ^ ^ ^ ^ Jpmi'^p+^n 



^p—^m—. 



a pm^O 

The lemma follows from H32() . □ 

Note that if e' < 1/4 then the last term in the estimate in Lemma iT.fil can be controlled by the 
positive term in Lemma 17.41 



8 Analyzing the quadratic Hamiltonian 

In this section we consider the main part of the Hamiltonian ^, namely the "quadratic" Hamil- 
tonian considered by Foldy. As in the previous section we have fixed z = 0, i.e., we consider the 
cube [— ^/2,£/2]^ and we omit z. The "quadratic" Hamiltonian consists of the kinetic energy K, 
from ((T^ and all the two-body terms with the coefficients w^^qq, ^oo,pq ^pO,Oq' ^Op,ij0^i^^ 
p,qj^O, i.e., 

N 

-f^Foldy = ~2~^' ~'~ ^ 5-/ ^P5,00^^ (^'^(^pe(^0e'(^qe''^0e + '^pe'^qe' %e' %e + ^Oe^Oe' ^pe'^qej ■ (^5) 

i=l PQT^O 

e, e' — lb 

In order to compute all the bounds we found it necessary to include the first term in (|34j) into the 
"quadratic" Hamiltonian. We therefore define 

Hq = Fpoidy + {^^ - ^ w^po.oo (ap+ao+ " a*_ao- + "*a-\-<^p+ - ao-«p-) • (36) 

p^O 

Our goal is to give a lower bound on the ground state energy of the Hamiltonian Hq. 
For any /c G denote Xk{^) = e*'^^x(x) and define the operators 

hl^ = {i''u^)-^l\l{VXk)a^± and h,,^ = {t'v^)-^'^ a^{VXk)al±, (37) 

where as before V denotes the projection orthogonal to constants. (On the subspace where = Q 
we set = and likewise when = we set = 0.) These operators satisfy the commutation 
relations (note that a^j^ commutes with aj^iVXk)) 



[hk+,hl'-] = n+,hl'-] = [bUK'-] = Q 



(38) 
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and^ 

In particular, we get 

K±X±] < ii'i^^)-'a*,^a,^ {VXk.VXk) < 1- (39) 

We shall now estimate the "quadratic" Hamiltonian below by an operator that is really quadratic 
in the operators and bf,^. 

8.1 LEMMA (The quadratic Hamiltonian). 

The quadratic Hamiltonian satisfies the lower bound 

Hq> hQ{k) dk - y2 w ooiop^a^ + a*_a ), (40) 

where 

hQ{k) = (2vr)-^ I''' j {bl^b,^ + b'Lk+b,k+ + K-bk- + bU-b-k-) 



+ ee'V^{bl,bk,,+bUeb-ke'+bleb-ke'+heb-ke')\ (41) 

e,e'=± 

Proof. We consider first the kinetic energy 

AT 

^^/Ci = ^^(a*+ag+ + ap_a^_)(Mp,/Cng) 

/I ^ I ^ 
|fc|2 _^ (^t6)-2 (ap+«q+ + ap-ag-)('Up,XA:)(Xfc,Ug)(iA; 

-3 7^ [ \^ 

2 j |fc|2 + (^t6)- 



(27r)-3i£ii / ; :„^^,^^ {a\{VXk)a+{VXk) + a*_{VXk)a^{VXk))dk. 



We get the first line in the expression for hq^k) if we use that 

aWPXk)a±i'PXk) > {v^)-^al{VXk)ao±al^a^{VXk) = e'blJ)k±- 

Since 



Wr,R{x,y) = {2tt) / Vr,R{k)Xk{x)Xk{y)dk, 

we have that 

-"^pg.oo = (2vr^)"^ / V;,/j(fc)(iip,PXfc)(^Xfc,"g)'iA; = (27r^)~3 / Vr,R{k){up,VXk){ug,VX-k) dk. 



If we use that K-,_r(^) = K-,r(— fc) and x(fc) = X{—k) we arrive at the expression for hq. The last 
sum in (IIU)) comes from commuting aQj^a^j^ to aQ_|_ao±. □ 
^The corresponding equation (25) in (LSo) has a typographical error. 
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8.2 THEOREM (Particular case of Bogolubov's method). 

For arbitrary constants A,B+,B- > and k £ C we have for each fixed /c G the following 
inequality involving the operators bj_f^j_ defined in 

A{bl+bf,_^ + + bl_b^_ + b*_^_b_^_) 

+K [Bfbl^ - b'J\_ + B'fb_,^ - B'Jh_, 

+K [b!/\^ - b'JX^ + B'/hu+ - B'J'bi, 

—ke' 



e,e'=± 



> -{A + B+ + B.) + ^{A + B+ + B-f - {B+ + B.)^ - \k\\ (42) 
Proof. Let us introduce 

4 = (^+ + i3_)-i/2(fiy'6i,+ - b'J%,_), 4 = iB++B.)-'/\B]/%,^+B'/%,^). 

These operators are analogous to the and df, discussed at the end of Sect. |21 According to 
and these operators satisfy 

K,di] = o, [c;,ci] = o, [d+,di]<i, [d_,d*_]<i, [c+,c;]<i, [c_,c*_]<i. 

The right side of (|42() can be rewritten as 

A{d*^d^ + d*_d_ + + c*_c_) 

+{B+ + s_) + d*_d_ + did*_ + d^d_ + {B+ + ^_)-i/2 (^^^^ ^ ^ ^(^^ ^ 



We may now complete the squares to write this as 

A{c\c^ +c*_c_) + D{d\ + ad^ + a){d\ + ad_ + a)* + D{d*_ + ad_^ + a){d*_ + ad^ +a) 
-Da^{[d^,dl] + [d_,d*_]) -2D\a\^ 

if 

D{l + a^) =A + B+ + B^, 2Da = B+ + B-, aD{l + a) = k{B+ + B-)^/^ . 



We choose the solution « = 1 + - ^(l + - L Hence 

Da^ = {B+ + B-)a/2 = \[a + B+ + B^- ^/{A + B+ + B.f - {B+ + ^_)2) , 

' ' Z)(l + a2 + 2a) " ^ + 2e+ + 2B_ - 2 1"! • 



□ 
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Usually when applying Bogolubov's method the operators b and b* satisfy the canonical com- 
mutation relations. In this case the lower bound given above is indeed the lowest eigenvalue of the 
quadratic operator. In our case the operators b and b* only satisfy the commutation inequalities 

We can now apply Bogolubov's method as formulated in the theorem above to the quadratic 
Hamiltonian Hq. 

8.3 LEMMA (Lower bound on the quadratic Hamiltonian). 

The quadratic Hamiltonian satisfies the lower bound 

Hq > -7m /'/o-^/^r^/^ - i (.+ - .-)^ ^)oo,oo " Ct'^.r \ (43) 
where Iq was defined in |3|). Likewise, 

^Foidy > -%,l/\u'/'r'/' - ct-'.r\ (44) 

Moreover, 



N 

Hpoidy > E " - Ct-'iyr\ (45) 



i=l 

Proof. We shall use (001). Note first that 

by ^ and the fact that R<1. 

By (grj) and Theorem lO with B± = v^, k = (z^+ - i^~)xik)i~^, and 



2(2^)3 / ' ' 4 |A:|2 + (£t6)-2' 



we have 



hQ{k) > (-(^ + r.) + + .^)2 - Z.2) _ (^+ _ ^-)2 |x(fe)|' 



i(27r)-3 ^5(A:) + /(fc) - (^(^(fe) + /(A;))' - <7(A;) 



l/2^ 



2(27r)3^6 



where 



Since 



2 |/fc|2 + (£t6)-2- 

g{k) = Attu {{k^ + R^)-^ - {k^ + r^y^) < 4Tru\k\-^ 
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and the expression g + f — {{g + f)^ — g^)^^^ is increasing in g for fixed / we have 

hQ{k) > -^(27r)-3 |^4Hfcr' + 7(A;)- ((4HA;r' + 7(A;))'-(47rz.)2|A;|-4y^'j 

If we now replace k by i'^/'^£~^/^^^ t^^k, in the ^-integration, and observe that 

= i^"^ J j Xi{x)Vr,R{x - y)Xe{y) dx dy = ^u^oo.oo^ 

we arrive at 

i(27r)-3 I 47r|A;|-2 + f{k) - ((47r|A;|-2 + f{k)f - {A'Kf\k\-^Y''^ dk (46) 



where 



1 = 



with 



' 2,,,. . „l/2^,,,^_i/2,_i2 



We now study the /-integral in more details. Since the integrand is monotone decreasing in 
/ we get an upper bound to / if we replace / by a lower bound. Let a = ^7^{^(i/^)~^/^t~^^. If 
1^1 < (8a)^/^ we shall simply replace / by 0. If |A;| > (8a)^/^ we replace / by the lower bound 

i|fc|2-a. 

Note that 

((47r|A:|-2 + l\kf - af - {A'Kf\k\-^f'^ = {\\kf + A-k + - (|A:|2 + ^',T\k\-'')aY'^ 
> {{\\k\' + A'.){l-m-'a)+a''f" 

If \k\ > (8a)^/2 then 4:\k\~'^a > a^(|fe|^/4 + 47r)"^ and hence in this case 
{Mkr + l\kf-ar-{4nf\kr)'/' 

> + ^-f' - m-'a CM + ^-f' + 2i\k\vf+A.y/^ - + 

> (i|A;|^ + 47r)^/^-a-Ca|A;r^-C|A;rV > {l\kf + Att)^^^ - a - Ca\k\-\ 
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Thus 

/ < ^(2vr)-^( J A7r\k\-Uk+ J A7r\k\-^ + ^\k\^ - {l\k\^ + 4tt)'^^^ + C\k\-^adk^ 

|fc|<(8a)i/2 |fc|>(8a)i/2 

< ^(27r)-3 / 4vr|A:r2 + i|fc|2- (i|A:|^ + 47r)^/^ dfe + Ca^/^ 

= io + c{ui)-^/H-' <io + Cu~'>/^f/H-^ur\ 

This gives the stated lower bound on Hq. The proof of ()44() is similar except that one should 
simply use Theorem 18 . 21 with k = 0. The proof of 1)45^ is the same as the proof of 1)44^ except that 
we, when proving the lower bound on -ffpoidy) replace IC by /C/2 and simply keep the other half /C/2 
in the lower bound. Note that this will of course change the constants in the lower bound. □ 

9 A-priori bounds on the kinetic energy, non-neutrality, and ex- 
citations 

In this section we shall estimate the kinetic energy (X^^Li^i); the non-neutrality — n~)^), 
and the excitations (z^ — n) in a state for which the energy expectation (H) is low. (Recall that H 
was defined in (|22|) and that we are omitting the subscript z.) More precisely, we shall introduce 
cutoffs 

R = ujit)-^e and r = £3/2(j,£)-i/2_ (47) 

and consider a state such that {H^J^) < 0. For technical reasons this assumption will be more 
appropriate later (see Lemma riU.l() . The Lemmas 16.317.117.217.417.51 fEa. and 17.61 together 

with (|45() in Lemma (8 . 31 control all terms in the Hamiltonian H^/^ from (|27() . If we combine all these 
estimates and assume that we have further cutoffs < r < r' < R' < R we obtain the following 

N N 
i=l i=l 

+i{(}oo,oo [(^"^ - - n+ - n-] - 4TTR''^i-^\n'^ -n-\{u-n) 

+^ E( E Jp>*p+%+ - - ^'"'(^ - ^) - 2^'^oo,oo ((^+ - ^-f + 1) 

-e'-'8nr^R'\n + l)(z. - n) - 2e' E( E " ^1-%-))^ ■ 

The prime on w'qq qq and J' refers to the fact that these quantities are calculated with the r' and 
R' cutoffs. If we choose s' = 1/8 and use '^00,00 — '^t^R''^^~^ we arrive at 

N N 
i=l i=l 



bound 

^r,R ^ 
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-^i^R''^ - Cz^2 (£-3/2^/1/2 ^ ^-3^/2) _ _ ^) _ cr^R'\n + l)iu-n + l) 

+Ho,oo(^+-^")'- (48) 

It turns out that the quantity I'i is important and that the main contribution to our energy 
estimate will come from boxes where is large (but not too large). We shall need the control on 
the kinetic energy, non-neutrality, and excitations only for boxes where v£ > ew(t)^. 

9.1 LEMMA (Bounds on non- neutrality and excitations). 

Let R and r he as in \4'^ - There is a constant Ci > such that if Ci£uj{t)'^ < ui and CiiV£^ < 
then for any state such that {H'^J^) < we have 

{^^') - Ce-^lH-''v''lH-^l\vlfl^ (49) 
l,v -n) <Ce~^IH-^{vifl^ (50) 

and 

((n+ _ n-f) < ^^-3/2^-2^(^^)3/2_ (5;^) 

Proof. We introduce the cutoffs 

R' = ae'/\,^i)-'/'i and r' = e-'/H^/\i^er'/\ 

where < a shall be specified below. We observe first that r < r' < R' /2 < R/2, which will, 
in particular, imply that Wqqqq > C~^R'^£~^. That r < r' follows from the assumption that 
ul > CieLo{t)'^ > Cie if we choose Ci > 1. That R' < R follows from the assumption Ciew(t)^ < ul 
if we choose Ci > a^. Finally, r' < R' /2 follows from the assumption that CiNi"^ < if Ci is chosen 
large enough depending on a. In fact, we get that R'/r' = ae^l '^{yt)^^l t'^l'^ > ae^^^N^^/'^i^^/^ 
since u < N . 

We see that (PU]) implies 

AT 



i=l 

From the estimate ()48() we hence obtain 

N 



1=1 

_C^.5/4^-3/4(^l + a-le-V2(^^)l/4 + ^-1/2(^^)1/4(1 + £-9/8^3/4(^^)3/8) ^ «2^(^^)-5/4^ 

-ct-^ue-\ 

We have here used that n < v and that we may estimate + 1 < 21^ since we assume that v >1. 
Since e~^!t^{vl) < e^^Ni^ < C^^ we may rewrite this estimate as 

N 

K'i > Y.^-flC.+ {he,t^"-Ci'^'--"c)er\u-n) + lw',,,,^^^^ 

i=l 
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-c{t-'^vr' + t-^H-^). 

If we choose a and Ci appropriately we see that the second term on the right side above is bounded 
below by C~'^el^'^{v — n). From this and the assumption Cie < Cet~^ < (recall that u){t) = 
Ct~'^) we easily obtain (jl^l) and We obtain ((SIJ) if we use that ■^00,00 ^ C~^R''^i'^ = 



10 Localization of the number of excited particles 

Although the bound on the expectation value (z/ — n) given in Lemma 19.11 is sufficient for our 
purposes we still need to know that {{u — n)^) ~ (1/ — n)^. We shall however not prove this for a 
general state with negative energy. Instead we shall show that we may change the ground state, 
without changing its energy expectation significantly, in such a way that the only occurring values 
of — n are bounded by the quantity on the right side of (|^U]) . 

To do this we shall use the method of localizing large matrices in Lemma lA.ll of Appendix 1X1 
Let ^ be any normalized wave function with v particles in the box {zf\ + [— ^/2,£, 2]^. We may 
write ^ = X]m=o '^m^m, where for all m = 0, 1, . . . , z^, ^'m, is a normalized eigenfunctions of n with 
eigenvalue m. We may now consider the + 1) x + Hermitean matrix A with matrix elements 

We shall use Lemma lA. II for this matrix and the vector if) = (cq, . . . , Cu)- We shall choose M in 
Lemma IA. II to be given by 

M = Integer part of £~^/'^r'^{ulf/'^ . (52) 

With the assumption in Lemma l9 . 1 1 that v£ > Ciew(t)^ we may assume (Ci so large) that M > 2. 
In particular, we may assume that M > C~^e^'^^'^t^'^{ul)^^'^ . With the notation in Lemma lA. II we 
have A = {ip,ATp) = (^', ff^'^^'). Note also that because of the structure of H^^^ we have, again 
with the notation in Lemma lA.ll that dk = i( k > 3. Since M > 2 this means that the second 
sum in (|71|) vanishes. We conclude from Lemma lA. II that there exists a normalized wave function 
^ with the property that the only occurring values of v — n belong to an interval of length M and 
such that 

(vi/,^;/,v) > {^,k!r^) - cM-\\d,m + \d2m)- 

We shall discuss di = di{'^) and d2 = d2{^) in detail below, but first we give the result on the 
localization of n that we shall use. 

10.1 LEMMA (Localization of n). 

Let R and r be as in J^Tj ) and M as in \5i3\) . Let Ci > be the constant in Lemma \y.l\ and assume 
that Ci£ < Cieuj{t)'^ < vl and CiNl^ < e^. Then for any normalized wave function ^ such that 

(m/,^;/,^) < -M-\\dim + \d2m) (53) 
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there exists a normalized wave function ^, which is a linear combination of eigenf unctions of v — n 
with eigenvalues less than CM only, such that 

[^MiR"^) > -M-2(|diW| + |d2W|). (54) 

Here di^'if) and (i2(^) are given as explained in Lemma \A.l\ 

Proof. We choose ^ as explained above. Then H54|) holds. We also know that the possible n values 
of ^ range in an interval of length M. We do not know however, where this interval is located. 
The assumption H53|) will allow us to say more about the location of the interval. 

In fact, it follows from (El and dSH) that {^^Hl'^^^ < 0. It is then a consequence of 

Lemma l9. II that {v — n)^'^ < CM. This of course establishes that the allowed values of — n 
are less than (C + 1)M (which we of course just write as CM). □ 

Our final task in this section is to bound di{^) and d2(^). We have that di{'^) = , H^/^{l)'i?), 
where H^^^{1) is the part of the Hamiltonian H^^^ containing all the terms with the coefficients 
ivpq^^jy for which precisely one or three indices are 0. These are the terms bounded in Lemmas 17.51 
and 17.61 These lemmas are stated as one-sided bounds. It is clear from the proof that they could 
have been stated as two sided bounds. Alternatively we may observe that H^^^{1) is transformed to 
—H^^^{1) by the unitary transform which maps all operators a*^ and a^^ with p / to —a*^ and 
—cbp^. This unitary transform leaves the estimates in Lemmas 17.51 and 17.61 invariant. We therefore 
immediately get the following bound on di{'^). 

10.2 LEMMA (Control of di{^)). 

With the assumptions in Lemma \lU.l\ we have for all e' > 

Likewise, we have that d2{^) = , H^J ^{2)'^) , where H^^^[2) is the part of the Hamiltonian 
-^ri? containing all the terms with precisely two Oq^ or two Oq^. i.e., these are the terms in the 
Foldy Hamiltonian, which do not commute with n. 

10.3 LEMMA (Control of da). 

With the assumptions in Lemma \1(J.1\ there exists a constant C > such that 

\d2{^)\ < Ce'^IH-'^v^I^I-^I\vI)^I^ + ^T^r'^R^ (*, {v - n)(n + 1)*) . 



Proof. We consider the unitary transform that replaces all the operators a*^ and Opj. with p ^ 
by —iOp^ and iOp^ respectively. Under this transformation the Foldy Hamiltonian H35|) changes 
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into an operator that differs from the Foldy Hamiltonian only by a change of sign on the part that 

N 



we denoted H'Jt{2). Since both operators satisfy the bound in (|44() we conclude that 



2=1 p<]¥=o 



If we use the representation ()31() we find using a Cauchy-Schwarz inequality that 



pq^O 



a \p^0 57^0 q^O p^O 

< ^ I ^ Jp0^p+'^0+ "^g0^g+'^0+) + ^ •^gO'^g-^O- JpO^-p-O-Q- 

X] ^^pg,oo(«p+(«o+ao+ + l)«g+ + ap-(aS-ao- + l)ag- 



The lemma now follows from (|32|) and from Lemma l9.ll □ 



11 Lower bound on the energy in a small cube 

Our goal in this section is to prove a lower bound on the expectation {H) = (4', H"^) of the operator 
H from (|22j) . We are aiming at the following result. 

11.1 THEOREM (Lower bound on H). There exists a constant Cq > such that for all 
< e < 1 and < t < 1/2 with CqNE^ < e^, and £~^t~'^i < Cq^ we have the estimate 

-Ki{e, t, N, £){nf/H-^/^ - {nf/H~^/^{N^/^i)~^/^ (c7eV6 + ^^^(e, t, N, , (55) 

where Ki and K2 are sums of a finite number of terms of the form Ce~°'t~^{N'^^^£YN~'^ with 
d>0. 

Note that the estimate is not formulated as an operator inequality. The lower bound is not a 
constant nor an expectation value, but a non- linear expression in (n). The reason that we give the 
bound in terms of an expression in (n) rather than in terms of v is that we shall later take into 
account the term T(S'o), i.e., the kinetic energy between boxes, which depends on (n) (see (|18p ). 

In order to arrive at an estimate expressed in (n) we shall make repeated use of the following 
result. 
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11.2 LEMMA. If {H) < we have the estimate 

(n) > v{l - C7e-i(iV2/5^)Ar-i/i5), 

In particular, we may assume that the constant Ci from Lemma \y.l\ is so large that if CiNf 
(i.e., (iV2/5£) < Cf^eiVViSj then 

V < C{n) and v^l^ < {nf/\l + Ce-\N^/^i)N-^/^^). 

Proof. This follows from Corollary 16.41 since e" we have used 

that < A^. □ 

We shall prove Theorem II 1 . 1 1 bv dividing into several cases. We first prove an a-priori bound 
on H. 

11.3 LEMMA (A-prori bound on H). We may assume that the constant Ci from Lemma \y.l\ 
is large enough so that if CiN£^ < e'^ and e~^t~^£ < C^^ there exists a constant C > such that 

H > -7-y"/oi^'/"r3/4 - Cu'/H-^/^ 1/4 + e-2(^£)3/4^^ _ c{t~'ur' + er^). (56) 

Proof. If = there is nothing to prove. We may therefore assume that u > 1. We now make the 
following choices for the cutoffs 

r = be'H^ and R = mm{ae^/^{ue)-^/^ ,uj{t)-^}e, 

where < a,b shall be specified below. We first note that r < ii if Ci is sufficiently large 
(depending on a and b). Indeed, either R/r = {a/b){e'^ij~^i~'^)^^'^ > {a/h)C\^'^ (since v < N) or 
R/r = b^'^uj{t)-^ee-'^ > C-^-H'^el-^ > b-^C-^Ci. 

We proceed as in the beginning of Sect. |3 but we now use Lemma 16.21 instead of Lemma 16.31 
and instead of (jlSI) • We then get that 
N 

H > ^-^(6-5)A,Neu-7e:yV^/'r3/^-C't-V-l-Cz.i?~l 
i=l 

-C7i/2(r3/2ri/2 + r^r^) + ^Wqo oo [(^^ " -n+- n'] - 'i-KRH'^\n+ - n-|(zy - n) 
'^^J2{Y1 '^pmK+"m+ ~ ap-«/.-)) - I^"^'^ - n) - 2e'woo,oo ((^^ - ^~)^ + l) 

If we Choose e' = 1/4, use that Wqq qq < 47ri2^£~^, and n < we obtain after inserting the choices 
of r and R that 

AT 

H > j;-^(.-5)A,Neu-%,(V^/'^-^/^-t^i-'^r'-^^-'a-iz.ri(z.£)V2 

i=l 
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We have here inserted the upper bound ae(z/£)~^/^£ for R. This is ahowed, since the only term 
which is not monotone increasing in R, i.e., ^i/R"^, can in fact, according to (|29j) . be ignored when 
R = u>{t)^^£. We have again replaced + 1 by 2i^. 
If we now choose 5 = e/2, and again use that 

N 

^ -^{e - (5)Aj,Neu > l'ye,te'^'^i~'^{i' - n), 

i=l 

we arrive at 

We obtain the results of the lemma if we make appropriate choices of a and b. □ 

We now use this estimate to control boxes with few or with many particles. We shall show that 
these boxes do not contribute to the leading order estimate of the energy. Thus the specific form 
of the first term in H56j) is not important. We may simply estimate it by Cv^^^i'^^'^. 

11.4 LEMMA (Boxes with few particles). 

Let Ci be the constant in Lemma \y.l\ and assume that vl < C\£u}{t)'^ , CiNl^ < , and e^^t^'^i < 
C^^. Then there is a constant C > such that 

H > -Ct-^^ur^ - Cel-^ - Ce-^t-^u. 

Proof. The result follows immediately from Lemma 111. 31 if we simply insert the assumed bound on 
u£ and use that uj{t) = Ct~'^. □ 



11.5 LEMMA (Boxes with many particles). Let Ci be the constant in Lemma \y.l\ Assume 
that vl > e-4(7V2/5^)io^ Q^^f < ^3^ £~H~^e < C{\ Then 

{H) > -^:}'^L^{nf/H-^/^-C{nf/H-^/^e~\N^'H)N-^'^'' 

-c{t-^^vr^ + t-^H-^). 

Proof. We may of course assume that {H) < and hence use the estimates in Lemma 111.21 We 
now use Lemma 111.31 and insert the estimates 

^5/4^-3/4(^^)1/4 ^ ^5/3^-l/3(^^)-l/6 < ^5/3^-1/3^2/3(^2/5^^-5/3 
^,5/4^-3/4(^^)3/4^ = z.5/3^-1/3^1/3^4/3 < ^5/3^-l/3(^2/5^)4/3^-3/15_ 

In the first inequality we used the assumption on v£. In the last inequality we simply used that 
iy< N. □ 
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We now restrict attention to boxes with eujit)'^ < < e^'^{N'^^^lY^ . In the next lemma we 
shall prove the lower bound on {H) under the restrictive assumption given in H57|) below. Finally, 
Theorem lll.ll is proved by considering the alternative case that H57|) fails. Let us note that, logically 
speaking, this could have been done in the reverse order. I.e., we could, instead, have begun with 
the case that (|57() fails. 

11.6 LEMMA (Lower bound on H: restricted version). Let R and r be given by i4'^ - Let 
Ci be the constant in Lemma\^ Assume that Cieuj{t)^ < i^i < e'^{N^/^iy° and CiNi'^ < e^. 
Then there exists a constant C2 > such that if 

£-^R^u{u - n) (57) 

\ a pn^O 



we have that 



h) > -^;}'ho{nf/H-'^/^-ct-^vr^ 

-C(n)5/^r3/4(e-i(Ar2/5£)Ar-i/i5 + ^-7/2(^2/5^)23/4^-3/ 
+e-9/2t-2(iv2/5£)8^r-i/5y 
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Proof. We may of course assume that {H) < and hence use Lemma 111.21 We may also assume 
that (H^^^) < otherwise we have from H29|) and the choice of r that 

{H) > -Ci.2(^-3/2^1/2^^-3^2)^_^^-3/2^5/4^-3/4(^^)l/2^3/4_^^^-l 
> _^^-3/2(^)5/4^-3/4(^-4(^2/5^)10)l/2(^2/5^)3/4^-3/10 _ ^^^-1 
= _^(^)5/4^-7/2^-3/4(^2/5^)23/4^-3/10 _ ^,^^-1^ 

which implies the stated estimate. It now follows from Lemma l9.1l that we have the estimate (|5U|). 

We again proceed as in the beginning of Sect.lUJ but we use (|44() instead of ()45() and ()29() instead 
of Lemma 16.31 or (|28|) (since now R = cj(t)^^^) and we choose 6 = e. We then get that 

+^-^00,00 [{^'^ - "~)^ - n+ - n~] - 47ri?^£~^|n+ - n~\{v - n) 

+^ Z]( X] - «p-«M-)) - - n) - 2e'tI}oo,oo ((^^ - n'f + l) 

If we use the assumption (|57j) and the facts that \n~^ — n~\ < u, n + 1 < lu, and t5oo,oo — 47ri?2^"3 
we see with appropriate choices of e' and C2 that 

{H) > -^-ji^i^y^lH-^i^ -C{t~^\i-^ ^t-'^H-'') 

-Ce-^l^v'lH^I^ - Cur^ - cu{tr^vr^ - C{u - nyH'\ 
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where we have inserted the choices of r and R. If we use ()5U|) we see that 

_ fi)j.l/2^-l < ^£-3/2^-2^5/4^-3/4(^^)3/4^1/2 < ^^-9/2^-2^5/4^-3/4(^2/5^)8^-1/5^ 

Hence we arrive at the bound in the lemma if we note, as above, that 

g-3/2j^7/4^1/2 < (^^5/4^-3/4^-7/2 (jy2/5£)23/4jy-3/10 

and use the result of Lemma 111.21 □ 

Proof of Theorem \ll.l\ We again let r and R be given by (|47|). According to Lemmas 111.41 and 
111.51 we may assume that Ci£uj{t)'^ < < £~^{N'^/^tj^^ . Then from Lemma lll.6l we may assume 
that 

r^R^v{v - n) (58) 

\ a p^jLO I 

We may also still assume that {H) < and {Hl'^) < (as in the proof of Lemma lll.6() . Hence we 
can use the results of Lemmas 19. II and 111.21 It is enough to prove an estimate of the type (|55|) for 
(H^^^) (the extra error terms were estimated in the first paragraph of the proof of Lemma lll.6p . 
We begin by bounding di and d2 using Lemmas I1U.2I and I1U.31 We have from Lemmas 19. II and 

nrni that 

1^21 < Ce~l/2^-2z.5/4£~3/4(^^)l/4 ^ cr'R'u{u - u) 

< Cv'^lH-^'l^ (e-^/H-^ue)^^ + t6e-3/2(^^)5/4^ < ^^6^-3/2^5/4^-3/4(^^)5/4^ 

where the last estimate follows since vi > ew(t)^. 

In order to bound di we shall use (|5H)l . Together with Lemma ri0.2l this gives (choosing e' = 1/2 
say) 

\di\ < Cr^R^u{{u - n) + 1) < C7t6e-3/2^5/4^-3/4(^^)5/4^ 

If the assumption (|53|) fails then 

{Hl[^) > -Ctl0e3/2^5/4^-3/4(^^)-7/4 > ,^^10^3/2^^-1 (^^(^^2 )-3/2 ^ _C7t22^^-l^ 

and we see then that the bound l|55|) holds. We may therefore assume that (|53|) holds. Thus from 
Lemma llO. II it follows that we can find a normalized wave function ^ , such that 

{v - n)$) < C7e-3/2^-2(zy£)3/2 and {v - nf^^ < Ce'^r^ivlf (59) 

and such that 

'M'K)>{^^Hl''^)-Ct''-r\ (60) 
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(Here as before {X) denotes the expectation of the operator X in the state 

In order to analyze (^^,H^^^^^ we again proceed as in the beginning of Sect.lHl but this time 
we use (|34p of Lemma 17.51 and (|43p of Lemma 18.31 We find that 

+ ^Wqq qq [(n+ - - n+ - n~] - iirR^i'^ln'^ - n~|(z/ - n) - ^u?oo,oo('^^ ~ ^~)'^ 
-e'-'87rr'R\n + -n)-2e'Y,{j2 Jp,i''p+%+ " ^1-%-))' ■ (61) 



We shaU use that for all e" > 



1 



and 

In'^ — n~\{v — n)< — in^ — n~S^ H iv — 



According to Lemma 17.41 

Z] Z] '^p/^('^p+^/i+ ~ '^p-'^^-)) — Z] '^pq,^iu{^p+'^q+^i>+^p+ ~^ '^p-'^q-^u-^^i- 
Oi p/J-J^O pq,fJ.UytO 



-2a*_|_a*_a^_a 



7.+) + - "^)'^ ^ 



The sum on the right is simply the second quantization of the two-body operator V ^Vwr^nV 
where V again denotes the projection onto the subspace of (([— £/2, ^/2]^) x {1, —1}) consisting 
of functions orthogonal to constants. Since Wr^R < this sum is bounded above by 



- n+f + {v- - n'f) r'^ <{v -n) 



(This can of course also be proved by directly estimating the sum). If we insert the above estimates 
into ()61jl and choose e' = l/e" we arrive at 

-CRH~^ {u + e"(n+ - n- f + (1 + pz) {{i^ - + l) + e"iy{iy - n)) 
-Cr-\^{u - nf + (1 + - n)). 

We now use (|59p . Lemma 19. II fwe may of course assume that (^', HI^r"^) < 0), and the choices ^ 
of R and r. We get 

-Cu{t)-H-^ {v + e"e-^/H-^u{utf'^ + (1 + jrr)e-h-\utf 
-Cl-^'\vlfl^ (^j,e-h-\vlf + (1 + ^)e-^/H-\ulfl^) . 
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Finally we make the choice arrive at 

If we now insert the bound vl < e~'^{N'^^^iy^ , use ()6Up and the result of Lemma lll.2l we arrive at 
a bound of the form for {H^^^) . □ 

12 The lattice approximation 

Given a map 5" : ^ M we define a function : — > M as follows. On any cube {^} + [— 1/2, 1/2]^ 
with /u E {(1/2, 1/2, 1/2)} + we set 

0(x)= ^ A,(x-;u)S(^ + r/2), (62) 

t| = ^ 

where 

A,(x) = (1/2 + TiXi)(l/2 + r2X2)(l/2 + T3X3). (63) 

Note that the 8 points /i + t/2 with |t| = \/3 are the corners of the cube {fi} + [—1/2, 1/2]^. The 
function cp is well-defined on and is continuous. Moreover, (j){cr) = S{a) for a G Z^. By a 
straightforward calculation we obtain that 

{Vcpf = T{S), (64) 

where we have defined 

^(^)= E ^(S(ai)-5(a2))2+ Y4^S{a^)-Sia2)f. 

kl-'^2l = V^ |CTi-o-2l = v^ 

Note that there is a constant C > such that for all maps S : — > R we have 



12.1 LEMMA. // Ai, . . . , A„^ > and YlY=i Aj = 1 we have for all Si, . . . , Sm > and all (3 > I 
that 

j=i j=i i=i i<i<i<'Ti j=i 

for some constant Cm,i3 > 0. 
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Proof. The first inequality follows from an application of Jensen's inequality. If we write Y = 
EJLi ^jSj then X]7=i Aj5^ - (ET=i ^J^jY = E7=i ^ji^j - Y'^)- The second inequality above 
then follows easily from \a^ — b^\ < f3\a — b\{a^~^ + b^~^), which holds for all a,b > 0, and 
\Si-Y\<ZT=i>^j\Si-Sj\. □ 

12.2 LEMMA. If (3 > 2 and if (p : ^ is the function constructed above corresponding to a 
non-negative map S : Z,^ ^ we have for all < 6 that 



[1-S)y S{af - C6-^^-^^T{Sf'^ < I ^^<y S{ 



and that 

S{af^ - 5T{S) -C6-\Y^S{afy\Y^S{af) < [ 0^/^ < J] S{af^ 
Proof. The functions Xr{x) in ()63() are non-negative and satisfy 

and that /[_i/2 1/2]^ \T{x)dx = 1/8 for all r G with |r| = \/3. We thus get from the inequality 
in Lemma nm that for ah /i G {(1/2, 1/2, 1/2)} + 



{/.}+[-l/2,l/2P 



> J E 5(/. + r/2)/^-C( E (5(A^ + n/2)-5(/. + r2/2))2)'^' J] + t/2)/3-i 



tGZ3 Ti,T2eZ'' 

If we sum these inequalities over /i (i.e., over cubes), use that each point in is the corner of 8 
different cubes, and that b^^"^ < b^)^^'^, when b^ > 0, we obtain the first inequality of the 
lemma. 

To arrive at the second inequality we instead use the Cauchy-Schwarz estimate 
E (5(/. + n/2)-5(^ + T2/2))2)'^' Y S{i, + T/2f-^< 

Ti,T2ez3 tGZ3 
|-ri| = \/3 = 



5 E {S{^^ + Tl/2)-S{^Ji + T2/2)f +C5-^ E ^(/^ + ^/^^ 

T]^,T2ez3 tSZ3 

|ti| = \/3 KNv^ 
|T2l = v/3 



2/3-2 
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For P = 5/2 we have 2(3 — 2 = 3. The second inequahty in the lemma then follows by summing 
over n as above, possibly redefining 5, and using that 



□ 



12.3 THEOREM (The Sobolev inequality for T). 

There exists a constant C > such that for all maps S" : — > M that vanish outside some finite 
set we have 

\S{aty^' <CT{S). 

Proof. Since T{S) > TdS"!) it is enough to consider only non-negative maps S. If we use the first 
inequality in Lemma 112.21 with /3 = 6 and 5 = 1/2 we obtain that 



^ \S{a)f <2 j (j){xfdx + CT{Sf 



on K" 

3 



Since : R'^ ^ M is a function of compact support we may use the standard Sobolev inequality 

xfdx <c(^j \V(t){x)\'^dx^ = CT{Sf. 

□ 

If S vanishes outside a finite set we get from the Sobolev inequality and the second inequality 
in Lemma 112.21 that 

^ S{af/^ - 5T{S) - C6-'(^Y1 '^(^)')' ^ /, "^'^^ ^ E ^(^)'^'- (^5) 



12.1 Application of the lattice approximation to the situation in Theorem 15.21 

We shall use the above results for the map 

So(a) = St {a) = r\^{n,) + l - 1). 

The expression T{Sq) is then exactly the term that appears in Theorem 15. 21 on localizing into small 
cubes. Note that Sq{(t) < {n^). Let i?!* : ^ M be the function corresponding to 5*0 as 

constructed in ()62() . Then from ()64() and Lemma 112.21 we get 

T{So) = j (V0)2 and f j (j)"^ < fS^iuf < ^ (n^) < N. (66) 
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Moreover, we get for all < 5 that 
CTGZ3 o-ez3 

< (1 _ 5)^7/4 f ^5/2 ^ ^ ^'5o(a)2)' + c6-^/^r'/\L/ef 



where we have used firstly that {n^) 7^ for at most (L/i)^ different points a and secondly the 
inequality 1)651) with 6 replaced by 6£~'^^'^. Likewise we get from the Sobolev inequality Theorem ll2.3l 

< CT{SQ){N'^'^Cf/^ + Cr'^''\L/lf. (68) 
Finally, note that if we define ^{x) = l-^l'^(\){xl I) we find that 

y (v$)2 = y (V0)2, y $5/2 ^ ^7/4 y ^5/2^ j^^=f jcf>^. (69) 

13 Completing the proof of the lower bound in Dyson's formula 
Theorem EH] 

In this final section we shall combine all the previous results to conclude the asymptotic lower 
bound in Theorem 12.11 From Theorems 14.21 111.11 and (|^7|) . 1)68(1 . and ((H^ we obtain, 

under the assumptions that CqNI^ < and e~^t~'^i < C^^ , that 

E{N) > inf (^A j (V$)^ -B j 5^/^^ - Z)iV^/^ 

where the infimum is over all functions > , such that J < N and 
A = {m,t/2 - 11%}'^ h5 - Ce^"" - CK^ie, t, N, i)) 

+i-6(^l/5^)-l^-l/5 ^ ^-1^-8^-2/5 

+77(7,-^/^/0 + E:i(e,t,iV,£))r3/2(iVl/5^)3(^2/5^)-15/4^-l/2 

+{e^'^ + K2{e, t, N, l)){N'^/^i)-^{N^I^LfN-^l^^ . 
We have used that = ^ and Y.a, u^^o < ^V^^- 
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We shall show that as ^ oo we can choose e,t,6 ^0 and N^/^£, N'^/^L ^ oo in such a way 
that A ^ 1/2, 5 ^ /o, -D ^ 0, e-^Ni^ 0, e'^t-^i 0, and i < L. 

Note that from ^ and (|21|) the limits of A and i? will follow if we can prove that Ki,K2 0. 
The error term D is of the form 

D = Ct-^{N'^/^L)-^ + Cr^N^/^iy^ + C{N^/^Lf{N'^'^l)-^ + K3 

where is a finite sum of terms of the form 

C6~''e-''t-^{N^^^e)'^{N^/^LyN-f, (70) 

where / > 0. Recall that Ki and K2 were sums of expressions of the same form. Note that 

e-3iV£3 = e-3(Ar2/5^)3^-i/5^ e^ir^^ = e-H-\N^/^e)N-^/^, and £/L = (A^2/5^)(^i/5^^-i^-i/5 

(which is required to be less than 1) are again all of the form ()7Up . 

Now first choose N'^/^i ^ 00 in such a way that (iV2/5^)d^-/ ^ Q foj. all ti^g 

occurring terms 

of the form Then choose N^'^L ^ 00 in such a way that C{N^/^Lf{N'^/^iy^ and 

{N'^l^lY{N^I^LYN-f for all the terms of the form Finally choose e,t, J ^ in such a 

way that all the terms of the form H70() and the first two terms in D still go to zero. Hence we have 
achieved all the limits as claimed. 

Dyson's formula now follows since by scaling we have that 

A j (V$)^ -B j ^ (2A)-3/5(5//q)4/5^7/5 j (y^)2 _ j ^5/2^ ^ 

where 

^x) = {2AIo/Bf/^N-^l^^^ [{2Ah/Bfl^N~^l^x^ 
satisfies / = iV"! / 5^ < i if J $2 < 

A Appendix: Localization of large matrices 

The following theorem allows us to reduce a big Hermitean matrix, to a smaller principal 
submatrix without changing the lowest eigenvalue very much. The proof can be found in |LSoj . 

A.l THEOREM (Localization of large matrices). Suppose that A is an N x N Hermitean 
matrix and let A^ , with k = 0,1, . . . , N — 1, denote the matrix consisting of the k^^ supra- and 
infra- diagonal of A. Let ip G be a normalized vector and set d^ = {ip,A^ip) and A = {ip,A'ip) = 
'Y2k=o ^k- (V' need not be an eigenvector of A.) 

Choose some positive integer M < N . Then, with M fixed, there is some n £ [0, — M] and 
some normalized vector (j) E with the property that = unless n + 1 < j < n + M (i.e., (p 
has length M ) and such that 

„ M-l N-l 

{cj,,Act>)<X + j^^k^\dk\+C^\dk\ , (71) 

k=l k=M 

where C > is a universal constant. (Note that the first sum starts with k = 1.) 
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B Appendix: Localization of the kinetic energy 

Our goal here is to prove a certain lower bound on the many body kinetic energy J^^^ — Aj acting 

N 

on the symmetric tensor product space 0gL^(M^). In order to state the bound we need to introduce 
some more notation. Let Xq denote the characteristic function of the cube [—i/2,£/2]^, in M^. Let 
Xz be the characteristic function of the cube {z£} + [—£/2,£/2]^, i.e., Xz{x) = Xq{x — £z). Let Vz 
denote the projection onto the subspace of L'^{{z£}+\—£/2, £/2]^) consisting of functions orthogonal 
to constants. We shall consider Vz as a projection in L^(M^). Let — A^^^^ be the Neumann Laplacian 
for the cube {z£} + [—£/2,£/2]^. Let aQ{z) be the creation operator 

a*,{z) = a*{£-^/^Xz), 

i.e., the operator creating the constant in the cube {z£} + [—£/2,£/2]^. Note that as before aQ{z) 

oo N 

acts in the Fock space 0^L^(M^). Products of the form aQ{z)aQ{z') or aQ{z')aQ{z) are however 

N=0 

N 

bounded operators on the space (S)sL'^{R^). Let for all z G the function Xz G C^{{zi} + 
{-£/2,e/2f) be such that \\d'^Xz\\oo < C(^t)-I"l and - ^^lloo < for some t with 

< t < 1 and all multi-indices a with \a\ < 3. We can now state the operator inequality we shall 
prove. Observe that that wc arc estimating the one-body kinetic energy operator (the Laplacian) 
from below by a many-body operator. 

B.l THEOREM (Kinetic energy localization). 

Let Q C M^. Let wi, . . . ,Wr G R^ and /3i, . . . , /3j. > 6e such that X^j.Xo = for all j = 1, . . . , r 
and such that for all v eM.^ we have that 

r 

i=i 

Then for all < s < t < 1 and all e > we have 

N 



(l + e + C{s/t)t'^ + Cs'/^ + Cs ) - A, 

3=1 i=l 

+ E f ( V«o(^ + ^3)%{^ + ^i) + 1/2 - \Jal{z)a^{z) + 1/2^ ] - ^ ^vol(O). 



Here all operators are considered in the sense of quadratic forms. 

Proof. By rescaling we may consider £ = 1. 

Note first that on the quadratic form domain of the Laplacian — A on we have 



/ (-ASu)^^ < / {-^in)dz = -A. 
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in the sense of quadratic forms. 

Let /s : R3 ^ R be the function fs{p) = p^/{p^ + s"^) for some s > 0. Then for all u G L2(R3) 
we have 

(27r)-^y" fs{p)\u{p)\^dp = J^x)(^-^^u) {x)dx. 



It follows that for all u G L^(R^) we have 



(27r)-6 I fs * \X\^(p)\u{p)\''dp = uix)X,{x) (^-^^{X,u)^ {x)dxdz. 



B.2 LEMMA. We have 

||/,-(27r)-V.*|^l'lloo<CsV2. 



Proof. We calculate 



2 ^2sin (ff3/2) y 



P3 / 

In particular Jj^i^^ \X{p)\^dp < Cr~^ for all r > 0. Since (27r)-3 / = 1, \\fs\\o<, < 1, and 

||V/s||oo < we have 

J\q\<r 

+ (27r)-3 / |/,(p)-/,(p-g)||X(g)|2dg<Csr + Cr-^ 

^|q|>r- 

The result follows from choosing r = s~^/^. □ 
We now write 

p^ = {27r)-Vfs*\X\\p)+p^[fM-i'^^r'fs*\X\\p))+p\l-fsip)) 
> (27r)-V/s * \X\\p) - Cs'/V+P\l - fsip)). 

For all u e L'^{R^) we have 

In other words we have the operator inequality 

- A > 1 A,dz + ^^f^ - CsV2(_A), (73) 

where the operator 

^ -A 
A, = Y,~d,X,-^^X,d, 

i=l 

is defined as giving the positive quadratic form 

^ r _A 

{u,Azu) = ^ / diu{x)Xz{x) _ _^ {Xzdiu){x)dx 
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with domain H^{{z} + [—1/2,1/2]^). With the similar notation the operator Yli=i~9i-^z9i is 
nothing but the Neumann Laplacian — A^^^. 

B.3 LEMMA. Given 9 G C°°{{z} + [-1/2, l/2]3) such that 9 is constant near the boundary of 
the cube {z} + [—1/2, 1/2]^. Assume moreover that ||5"6'||oo < Cjij"'"' for some t with s < t and 
all multi-indices a with \a\ < 3. Then 



[[A,,0],O] > -C(5/t)(-A£j - C{s/t)t-'X^ 



Proof. We calculate 



-A 



-A + s 



Xzdj. 9 



-2 



and 



-diX, 



A 



-A + s- 



rX^dj. 9 



-diX, 
-{di0) 

-2diX, 



-A 



-A + .s-2 

-A 
-A + s-2 



{di9) - diX, 

X,di 



-A 



-A + s-2' 



X,di 



-A 



-A + s-2' 
-A 



{di9) - 2{di9) 



A 



-A + s-2 



We also calculate 



-A 



-A + s-2' 



-A + s- 



{di9) - diX, 



[-A,9] 



-A + s-2^ 
A 



-A + s-2' 



X,di 



X,di. 



1 



-A + s-2 



and 



-A 



-A + s-2' 



-2 



1 



-A + s- 



+2s" 



1 



Hence 



Likewise 



-A 



-A + s-2' 



■A + s-2 

< Cs/t and 



[[-A,^],^] 
[-A,^]- 



1 



-A + s-2 

1 



-A + s- 



[9, -A] 



-A 



-A 



-A + s-2 



,di9 



,di9 



A + s-2' 
< C{s/t)H-\ 



-A + s-2 

< C{s/tf 



Using these estimates we find that 

-A 



-A + s-2' 
< -{t/s)diX, 



{di9) + {di9) 
-A + s-2'^ 



-A 



■A + s-2 
-A 



,0 



X,di 



-A + s-2' 



x,di + {s/t)idi9f 



< -C{s/t)diX^di + C{s/t)t~^X: 
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and 



-A + s 



3 1 

< E -9jidiO)—^—^idie)dj + C(s/t)r2x, 



3 



With = — we may write 



Here we have used that Xz~^zXz = Xz-z^p^Xz^ where —A is the Laplacian on M^. 
Clearly we have VzAzVz = Az- Thus 

Az > VzXz }^fl_^ XzVz - C(g/t)(-Ai,l) - C{s/t)t-^Vz 
^ -P^Xz }~fl,, XzVz - Cis/t)t-\-A'^l), 

(z) 

where we have assumed that t < 1 and used that — AJ^^^ > CVz- 
We turn now to the term — As~^/(— A + s^^) in (f7H|). 

B.4 LEMMA. For all u G we have that 

Proof. Write 6 = 2 f Pj j ■ On one hand we may write 



-A + s- 



□ 
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On the other hand we have 



^ J Pj \u{x + Wj) - u{x)\'^ dx = ^(27r) ^ J Pj\u{p)\'^\exp{i{p,Wj)/2) - exp{-i{p,Wj)/2)\'^dp 

= J^(2vr)-3y" P,\u{p)\Hsm\{p,w,)/2)dp. 



We thus simply have to prove that 



,2^-2 



p s 



p2 _j_ g 2 



+ bsp"^ > 4/3j sin2((p, tt;j)/2). 



If 2|p| < 6^/2^-1/2 j^g^yg 



p2 _|_ g~2 



> > (1 - bs)p^ > ^ 4/?,- sm2((p, Wj)/2) - hsp\ 



If \p\ > 6I/25-1/2 

we have 



bsp^ >b^>Y^ 4/3j sin2((p, u;j)/2). 



□ 



Since for ah y G M , (Xy+o-)o-£23 is an orthonormal family in L (M'^), we get from Parseval's 
identity that 

/ \u{x + Wj) — u{x)\'^ dx > / Xy^a{x){u{x + Wj) — u{x))dx . 

If we integrate this over y £ [0, 1]^ we obtain 

J \u{x + Wj) — u{x)\'^ dz > J J Xz{x) {u{x + Wj) — u{x)) dx dz. 

The quadratic form on the right of this inequality corresponds to the operator whose second quan- 
tization has the form 

j {a*Q{z + Wj) - a*Q{z)) {aoiz + wj) - aQ{z))dz. 
It is then clear that the main inequality (|72|) follows from the following lemma. 
B.5 LEMMA. Let a\ and he two commuting creation operators. Then 



{al - a*)(a2 - a^) > J a^a^ + 1/2 - J a\a^ + 1/2 - 1. 



Proof. We calculate 



(02 - a*)(a2 - Oi) = 0202 + aloi - {ala^ + 0^02) 
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We have 

= 2a2aia\a2 + 2a]|a2a2'^i = 40202*^101 + 2aja]^ + 20202 

= 4(0^02 + l/2)(o^ai + 1/2) - 1 < 4(0^02 + l/2)(a^^ai + 1/2). 

Since the square root is an operator monotone function we have that 



a^Qi + 0^02 < 2 J O2O2 + 1/2 J ala^ + 1/2 



(02 - Oi)(o2 - > 0302 + a5;oi - 2^/0302 + 1/2-^/ (o^o^ + 1/2 

2 



Hence 



0^02 + 1/2 - Y'a^ai + I/2J -1. 

□ 
□ 

We shall use Theorem IB. II with wi,W2, . . . ,Wr being the vectors a G Z'^ satisfying that jcrp = 2 
or Icrp = 3. Note that there is a total of 20 such vectors. We define 



Pa 



12 ' I I — 



We then have for all G R'^ that 



\a\ = V2 or \it\ = V3 1ct| = v^ \tT\^V3 
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